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Abstract 

We present quasi-analytical and numerical calculations of Gaussian wave packet solutions of 
the Schrodinger equation for two-dimensional infinite well and quantum billiard problems with 
equilateral triangle, square, and circular footprints. These cases correspond to N = 3, N = 4, 
and N —* oo regular polygonal billiards and infinite wells, respectively. In each case the energy 
eigenvalues and wavefunctions are given in terms of familiar special functions. For the first two 
systems, we obtain closed form expressions for the expansion coefficients for localized Gaussian 
wavepackets in terms of the eigenstates of the particular geometry. For the circular case, we discuss 
numerical approaches. We use these results to discuss the short-time, quasi-classical evolution in 
these geometries and the structure of wave packet revivals. We also show how related half- well 
problems can be easily solved in each of the three cases. 
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I. INTRODUCTION 



The use of wave packets to analyze the dynamics of quantum mechanical systems is 
an increasingly important aspect of the study of the classical-quantum interface.— Popular 
software packages^ can help students visualize the evolution of quantum states (in contrast 
to the more typical time-independent stationary state solutions seen in most textbooks) 
by allowing students to change parameters (such as the initial width of a wave packet). 
Such visualizations are potentially important because recent studies,— ^ have suggested that 
student understanding of time-development in quantum mechanics is rather limited. 

Several authors have discussed pedagogical descriptions of various one- dimensional (ID) 
quantum mechanical problems using a wave packet approach, studying such topics as trans- 
mission and reflection from square barriers , 5i^7i? linear potential steps,^ "bounces" off infinite 
walls , 10 ! 11 and bound state wave packets in single square wells in either position spaced or 
momentum space,— in double wells , 14 1 15 and in systems of relevance to solid state physics . 16 ' 17 
More recently, wave packet revivals have been studied in the context of familiar ID quan- 
tum mechanical bound state systems such as ID infinite wells and the "quantum bouncer" 
problem 19 - 20 - 21 - 22 - 23 - 24 - 25 - 26 - 27 - 28 

An obvious extension would be the study of wave packet propagation in two-dimensional 
infinite wells of various shapes (quantum billiards). In this paper we will discuss the construc- 
tion of Gaussian-like wave packet solutions in three simple two-dimensional (2D) geometries, 
namely, the square, equilateral triangle, and circular infinite wells. We note that measure- 
ments of conductance fluctuations in ballistic microstructures^ 9 - have been tentatively used 
to identify features in the power spectrum with particular closed orbits in a circular and 
stadium billiard. More recently, the realization of billiards^ - with ultracold atoms in arbi- 
trarily shaped 2D boundaries confined by optical dipole potentials has allowed the study of 
various chaotic and integrable shapes such as the stadium, ellipse, and circle. 

The study of the energy eigenvalues and eigenfunctions in general 2D billiard systems 
can be used to probe quantum chaos, but in the cases we discuss here, the corresponding 
integrable classical motions are easily obtained and can be compared to, and contrasted 
against, the short-time quantum development. In each case, we make use of closed form 
solutions, and for the first two systems we show how to obtain simple expressions for the 
expansion coefficients of a gaussian wavepacket in terms of the eigenstates of the particular 
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system; for the circular case we illustrate numerical results. We then use our results to 
discuss the short-time propagation of wave packets, compare them to classical orbits, and 
investigate the structure of longer-time quantum revivals. 

The infinite square well provides an easy introduction to the subject, and the cases of 
the circular— and equilateral triangle^ quantum billiards have recently been examined in 
the research literature. We will make appropriate use of the results in Refs. and ^2, but 
our discussions will focus on more pedagogical aspects, including other special cases of 2D 
billiards obtained by "foldings" of these cases to make connections to the degeneracy of the 
energy eigenvalues and the symmetry structure of energy eigenstates in simple geometries. 
The comparison of these three exactly soluble 2D infinite well problems and their quantum 
dynamics is a major thrust of this paper. Note that these systems correspond to N = 4, 
N = 3, and N — > oo regular polygonal billiard footprints, and it is an interesting open 
question how the exact quantum revivals we find for the N = 3 and 4 cases change to the 
approximate ones for the N — > oo (circular) case. 

We begin in Sec. II by reviewing the basic properties of Gaussian wave packets and the 
general time-dependence of ID bound state systems, before discussing the ID infinite well in 
Sec. III. In Sec. IV we consider the general time-dependence of systems with two quantum 
numbers, as well as the square billiard. We will also be able to use this geometry to examine 
the special case of the 45°-45°-90° triangular billiard obtained from the square well by 
folding along a diagonal. In Sec. V we briefly review the properties of the quantum solutions 
to the Schrodinger equation in an equilateral triangle (60°-60°-60°) infinite well (one that is 
not frequently discussed in the pedagogical literature) and discuss both the short-time quasi- 
classical propagation of wave packets as well as derive the exact quantum revivals that are 
present in this system. We also note that half-well solutions for the 30°-60°-90° triangle 
are easily obtained by folding along an axis of symmetry. Finally, we turn our attention in 
Sec. VI to the case of the circular billiard and again discuss how information on both the 
classical periodicities and possible quantum revivals is contained in the energy eigenvalue 
spectrum. In Sec. VII we present conclusions and some avenues for future investigation. 

Although this work will focus mostly on the formalism required to construct Gaus- 
sian wave packets in these three geometries in the most efficient ways and will be re- 
stricted to static images of position-space probability densities, we have also produced an- 
imations of the short- and long-time development of such wave packet solutions in both 
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the ID and 2D infinite wells, illustrating the short-time, quasi-classical propagation and 
the structure of quantum revivals (and fractional revivals). These animations are avail- 
able at <http://www.phys.psu.edu/~rick/QM/qm.html> and include visualizations of the 



position- and momentum space probability densities, quantum expectation values compared 
to classical trajectories, and plots of the autocorrelation function for several different initial 
conditions for ID and 2D Gaussian wave packets. 

II. ONE-DIMENSIONAL BACKGROUND 
A. Gaussian wave packets 

For each of the geometries we consider, we will use initial Gaussian wave packets cor- 
responding to standard free particle solutions of the Schrodinger equation. Although such 
solutions do not explicitly satisfy the boundary conditions at the wall(s), as long as they are 
sufficiently far from any infinite wall boundary, the resulting error can be made exponen- 
tially small, and the solutions can be reliably expanded in eigenstates of the billiard system. 
The initial wave forms we will use in one dimension are given by 

(p\G) = <f> G (p,0) = ^l e - 2 (p-po)V2 e - ipxo /n ^ (1) 

which corresponds to 

(x\G) = Mx, 0) = _L_ e -^-^) 2 / 262 e^^-^)^ , (2) 

where b = ah. These two equivalent representations of the same Hilbert state vector, \G), 
are, of course, related by the Fourier transform 

i r+oo 

VfcfoO) = / 4> G (p,0)e ipx/h dp. (3) 



The initial expectation values for this general state are given by 

b 2 b 

(x) = x , (x 2 ) = x 2 + —, Ax = -j=, (4) 

and 
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We can easily vary the initial position xq, momentum po, and width Axo of the packet. The 
expectation value for the total (kinetic) energy is given by 

^) = hp 2 ) = Uvl + ^). (6) 



2^ ' 2/iV^ u 26V' 

We will use the notation /i for the particle mass to avoid confusion with various quantum 
numbers. 

If we expand such an initial state in terms of the energy eigenstates u n (x) of a general 
ID bound state potential, we have 

oo 

V>g(z, 0) = ^a n w„(x), (7a) 

n=l 

with 



u* n (x) ipc{.x } 0) dx. (7b) 
The resulting a n must satisfy two constraints, namely, 

oo 

= 1 (8a) 



n=l 

oo 



J2K\ 2 E n = (E) = ^(pl + ^), (8b) 



n=l 



where E n represents the quantized energies of the potential. Equation (jSJ) is very useful 
for numerical checks on the expansion coefficients. (We note that even if the subsequent 
time development of the quantum wave packet involves "collisions" with an infinite wall, 
if the wave packet is initially normalized to unit probability, it will remain so through its 
subsequent time development. Once we normalize the wave packet, unitarity ensures that 
we can "set it and forget it.") The resulting time-dependence is then given by 

oo 

i/)(x,t) =Y,anU n (x)e- lE ^ h , (9) 

n=l 

and we kno w 11 ! 20 that for short times (before the initially Gaussian packet has a chance 
to "collide" with any of the infinite wall boundaries), its behavior will resemble that of a 
free-particle solution of the form 

\^(x,t)\ 2 = ^— e -(—o-p t/^/b^ (1Q) 

bt\/n 
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where 



with 



b t = ah^l + (i/t ) 2 , (Ha) 
Ax t = Axo^/TTJtJW, (Hb) 



t .^=(|)A^ (12) 



which defines the time scale for wave packet spreading. 



B. General time-dependence 

For the expansion of a general wave packet as in Eq. Q , we typically expand the energy 
eigenvalues about the central value of the quantum number n (assumed to be an integer 
for simplicity) used in the construction of the wave packet, namely, 

E{n) « E{no) + E'{n ){n - n ) + -E"{n ){n - n ) 2 + \E'"{n ){n - n ) 3 + • • • (13) 

We can express the time-dependence of each quantum eigenstate as 

e -iE n t/n = exp (_ z /^[ jE ( no ) t + („ _ n )E'(n )t + \{n- n ) 2 E"(n )t + ...]) 
= exp(— iuj^t — 1ixi{n — rio)t/T c \ — 2ni(n — rio) 2 t/T Tev 

-2m(n-n ) 3 t/T supci + ...), (14) 

in terms of which the classical period, quantum mechanical revival, and superrevival times 
are given respectively by 

r - - W^l (15) 

Jrev \E»(n )\/2 
T = 2nh (17) 

The first term of Eq. (JT4"Jl is an common overall phase that gives the same (trivial) time- 
dependence to each state, while the second term describes the short-term, classical peri- 
odicity of the bound state system because it gives e -27r H n -"o) — \ f or eac h state after one 
classical period, T c \. (An alternative derivation of the form T c \ = 2irh/\E f (no)\ is given in 
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Appendix A using a WKB approach.) Bound state wave packets (with the exception of the 
harmonic oscillator potential for which the energy spectrum is exactly linear in n) will typ- 
ically spread significantly after a number of classical periods, entering a so-called collapsed 
phase, only to reform later in the form of a quantum revival, in which the spreading reverses 
itself (completely in the case of the infinite well, partially in the case of a general potential 
well) and the wave packet relocalizes. The revival time, T rev , describes the (typically longer) 
time scale for that behavior. 

A standard tool for analyzing the short-term, semi-classical periodicity, and the long-term 
revival structure is the autocorrelation function defined by*^ 

/+oo 
ip*{x,t)^{x, 0)dx 
-oo 

/+oo 
<P*(p,t)<p(p,0)dp (18) 
-oo 
oo 

= K\ 2 e lEnt/ \ 

n=l 

which measures the degree of overlap of the initial wave function with itself at later times 
(in both position- and momentum space). 



III. ONE-DIMENSIONAL INFINITE WELL 



For the ID infinite well of width a, the energy eigenvalues and eigenf unctions are given 

by 

E n = ^ = - ^E (19) 

u n {x) = y-sin(-— ) . (20) 

The classical period obtained from Eq. (|15|) (using the identification fiv n = p n = rink/ a) is 
given by 



dn fxa 2 a \ fi 
which implies that 



dE n h 2 7i 2 n h% ( p n \ hirv n 

" ' ( 21 ) 



Tt = ™ = *, (22) 
\dE/dn\ v n 



as expected. The revival time is given as 



4/ia 2 2nh 

= ~Tm~ = "aT ' (23) 



and the revivals are exact in this case because e~ tEnTlcv ^ h = e -2 ™ 2 = 1 for all values of n. 
Because of the quadratic n-dependence, there are no higher order time scales beyond the 
classical period, T c \, and the quantum revival time, T rev . 

The expansion coefficients for a general Gaussian wave packet of the form in Eq. (J2J) are 
given by 

a n = [u n (x)][ip G (x,0)]dx , (24) 
Jo 

where we can do the integrations numerically over the finite (0, a) interval. We show in 
Fig. 1 the results of such a numerical computation, both for the normalization, Y^m=i \ a n\ 2 , 
and the average energy, Yl^=i E n \an\ 2 , for zero momentum (po = 0) wave packets as we vary 
the initial position, x$: we use here (and elsewhere) the numerical values 

h = 2/i = a = 1 , (25) 

and show results for values of b = l/(2vT0) (1/vTo) corresponding to initial spatial widths 
of Axo = 0.05 (0.10) as shown by the solid (dashed) curves; we used the 40 lowest-lying 
energy eigenstates in the expansion. As expected, as we move the initial wave packet near 
the edge (and beyond), the probability of being inside the well decreases uniformly, and the 
wider wave packet "senses" the wall first, as shown at the bottom of Fig. 1. Initial wave 
packets that are sufficiently far from either wall can therefore be reliably "supported" inside 
the well, despite the fact that they do not explicitly satisfy the boundary conditions at x = 
and x = a. 

Perhaps more unexpectedly, the average energy of these wave packets grows rather dra- 
matically (from the Po = value of (E) = (h 2 /2b 2 )/2/j,) as the packet nears the wall. This 
behavior is due to the expansion in Eq. (J7|) attempting to reproduce the sharply discontin- 
uous jump enforced by the edge of the well at x = a, which requires the inclusion of many 
shorter wavelength and hence higher energy/momentum component states. In fact, for a 
sufficiently poorly behaved discontinuous function, the evaluation of the kinetic energy in 
such a way may be ill-defined,— and increase without bound as the number of component 
states is increased. 

To avoid all such difficulties, we will henceforth assume that the initial Gaussian is suffi- 
ciently contained within the billiard footprint so that we make an exponentially small error 
by neglecting any overlap with the region outside the ID or 2D well, and may ignore any 
discontinuity at the wall. In practice, this condition only requires the wave packet to be a 
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few times Axq, = b/y/2 away from an infinite wall boundary. We can then extend the inte- 
gration region from the finite (0, a) interval to the entire ID space, giving the (exponentially 
good) approximation for the expansion coefficients 

a n = I [u n {x)\ [ip G (x,0)] dx 

J in-KX /a^—b 2 (p +rnrh/a) 2 /2h 2 ^-in-KX /a^— b 2 (p — nnh/a) 2 /2h 2 j (26) 




because we can write 



sm = i(e m »/« _ e -mW«) j (27) 

and we can perform Gaussian integrals such as 

hoo [— 

e~ ax2 ~ bx dx= J-e b2/4a (28) 
V a 

in closed form. This expression is very useful because it can speed up numerical calculation 
involving the expansion coefficients such as the evaluation of the autocorrelation function in 
Eq. (fTHj). It also accurately encodes the sometimes delicate interplay between the oscillatory 
pieces of the Gaussian (e~ ipoX / h ) and the bound state ( e ±m7nr / a ) wavefunctions, which can be 
difficult to reproduce in a purely numerical evaluation, and it does so in a way that is valid 
for arbitrarily large values of po, where the integrand would be highly oscillatory. Finally, 
Eq. ()2(j)l nicely illustrates how ipc{x, 0) and the u n (x) must not only have an appropriate 
overlap in position space, but also must have an appropriate phase relationship between 
their oscillatory terms. This phase connection leads to the exp(— b 2 (po ± nnh/a) 2 /2h 2 ) 
terms, which can be understood from a complementary overlap in momentum space. (We 
discuss this point in more detail, as well as the evaluation of the d n in momentum space in 
Appendix B.) 

We can make immediate use of Eq. ()2fi)l by considering zero momentum (po = 0) wave 
packets; this case corresponds to placing an object "at rest" inside the infinite well potential. 
For such cases, the only natural periodicity in the problem is the revival time in Eq. (|23j). 
because there is no classical periodic motion. In this special case, the expression for d n in 
Eq. ()26|) simplifies even further to 



i^e- bW ^ sin( !^), (29) 
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which shows that for several special values of xq in the well, several of the expansion coef- 
ficients will vanish for obvious symmetry reasons. For example, for xo/a = 1/2, all of the 
even n = 2, 4, 6, . . . coefficients are zero and the only non-vanishing terms in the expansion 
are the odd ones (n — 2k + 1) which have energies of the form 

E k = ^(2k + l) 2 = E (4k 2 + 4k + l) = E + 8£ [^±12] . (30 ) 

The first term in Eq. (J30|) contributes only to the same overall phase to the time-dependence 
of each term. The second term is of the form 8Eq times an integer and leads to revival times 
that are 8 times shorter than the standard T rev = 2nh/E in Eq. (}2*3"|) . Similarly, for the 
cases of x^/a = 1/3, 2/3, the a n with n = 3k vanish, leading to special exact revivals at 
multiples of T rcv /3 for these two initial locations. (It is also possible to construct odd-parity 
wave packets that have different patterns of special revivals at other locations in the well.) 

To illustrate both the revivals and fractional revivals in special cases, we plot in Fig. 2 
the autocorrelation function, |t4(£)| 2 , for one revival time, T rev , for different initial central 
positions, Xo/a, of the zero- momentum packet from the center to as near one edge as we can. 
We note that in each case there is an exact revival at T rev = 27rh/E , while for the special 
cases of xo/a = 1/2 and 2/3, we find exact revivals at shorter time intervals as expected. 
For xo/a = 0.8, we also notice large partial (but not exact) revivals at 0.4T rev and 0.6T rev 
for similar reasons (because sin(4ri7r/5) vanishes for n = 5k.) 

We next show in Fig. 3 the effect of "turning on" momentum values for an i = a/2 
wave packet. For p = (top line), we have the special pattern of exact revivals at multiples 
of T rev /8 noted above, due to the vanishing of the even expansion coefficients (shown in 
the corresponding |a n | 2 versus n plot in the right column). For a small, non-zero value 
(po = 37r, second row), only the exact revival at T rev remains, because the even expansion 
coefficients are no longer forced to vanish. The autocorrelation function decreases somewhat 
more rapidly from its initial value than in the po = case, because the particle is slowly 
moving away from its initial position, in addition to spreading out. 

For still larger values of momentum, such as po = 40tt (third row), we see obvious evi- 
dence for the classical periodicity and the first appearance of fractional revivals^ at rational 
fraction multiples of T rev . The corresponding a n now exhibit a more obvious Gaussian shape, 
with a spread, An, which is related to the momentum variable by p n = mrh/a, so that 

Ap = — An, (31) 



2Ax y/2b a 
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which gives An = a/2irAxQ. (This relation can be useful in deciding how many values of a n 
to include in a numerical evaluation of quantities such as A(t).) For even larger momentum 
values (see the p = 4007T case at the bottom, for example), the classical period becomes 
much shorter than any obvious fractional revival time scale, and the shape of the expansion 
coefficient distribution is unchanged (same An), but simply shifted to higher values of n. 

These results, especially the expression for the expansion coefficients in Eq. (J26|) . can be 
used in the square billiard we consider in Sec. IIVI because that geometry is clearly separable 
into two ID problems. 



IV. SQUARE BILLIARD 



A. Time-dependence for two-dimensional systems 

Systems with two quantum number o 35 i 36 with energy E(ni,n 2 ) offer richer possibilities for 
both semi-classical periodicities as well as wave packet revivals. The extension of Eq. (fE?|) 
gives two possible classical periods, namely 

(m) _ 2nh 
ld ~ \dE{n x ,n 2 )/d ni \ [6 a) 



T (n 2 ) = 27Tk 

cl ^ \dE{n 1)n 2)/dn 2 y 



and the long-time revival structure typically depends on three possible times given by 



2nh 

(l/2)|d 2 E^n 2 )/dn\\ 



Tt ] = T, ~T7^ (33b) 



(l/2)\d 2 E(n 1 ,n 2 )/dn; 
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T (ni,n 2 ) _ ^ (on \ 

\d 2 E{n 1 ,n 2 )/dn 1 dn 2 \ ' 1 ' 

and the revival structure depends on the interplay between these three times. 

Classical closed or periodic orbits, with periods given by T^°\ are reproduced when the 
two classical periods are commensurate, namely, when 

pTj; i} = = rf 2) • (34) 



The relation ()34|) can be interpreted as arising from the beating of the two classical periods 
against each other.— 
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B. Time-dependence for the square billiard 



For the infinite square well (with dimensions L x x L y = a x a), the problem simplifies 
to two copies of a single ID infinite well because of the separability of the potential. For 
example, the energy eigenvalues, E(n x ,n y ), and eigenstates, W( nx ,n y )( x ,y), are given by 



E(n x ,n y ) = — — , (35a) 



and 



W{n x ,n y ){x,y) = U(n x ){x)u(n y ){y), (35b) 



where the appropriate quantum numbers. The revival times are given 

by Eqs. (j33aj) and (|33b|) . and are simply related to each other via 

rp(n x ) _ _ T (n y ) /oe\ 

rev »- -'rev ' \" u / 

with no cross-term present. Therefore, the quantum revival structure also is very simply 
related to that of the ID infinite well, including the possibilities of special symmetric revivals 
for zero-momentum wave packets at particular locations, such as (xcbZ/o) — (°/2,a/2) and 
(a/3,2a/3). For rectangular infinite wells with incommensurate sides (L x /L y ^ p/q), the 
structure of the revival times may be more complex . 35 ' 36 

Compared to the simple back-and-forth motions in the ID infinite well, the closed or 
periodic orbits here are more interesting and the nontrivial path lengths for closed orbits, 
L(p, q), in the square billiard can be readily deduced from simple geometric arguments , 37 ' 38 
and are given by 



L(p,q) = 2avV + q 2 (37) 

where 2p and 2q count the number of "hits" on the horizontal and vertical walls respectively, 
before returning to the starting point in phase space. The corresponding classical periods 
for such closed trajectories would be given by 

T cl =-^~> ( 38 ) 

where vq is the classical speed. Such orbits can be produced by point particles in the 2D 
billiard, starting from any initial location, (x , y ), inside the box, provided they are pointed 
appropriately, namely in the tan(#) = q/p direction. The values of 9 = tan _1 (g/p) and 
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^ci l T = a/p 2 + (where r = 2a /vq is the period for the simplest, back-and-forth closed 
trajectory) for many of the low-lying cases are tabulated in Table I. 

The condition for periodic orbits in Eq. (J3"4^ can be implemented very easily in this case, 
and we will examine its derivation in detail. For such closed orbits to occur we require that 



/ 2/ia 2 x (n ) 



cl 



qT (n y ) 



, 2/ia 2 , 
Kim,, 



(39) 



or n y = n x (q/p). If we substitute Eq. (p?9"j) into Eq. (J33j) . as well as equate the quantized 
total energy, E(n x ,n y ), with the classical kinetic energy, we obtain 



2^2 



E(n x ,n y ) 



n 2 7r 

fr 2 7T 2 

2/ia 2 
2/ia 2 



2 . 2/9V 

n 2 x (p 2 + q 2 ) 

p2 



(40) 



or 



71; 



V 



and 



fiavo 



HIT ^Jpl + g2 ' 



so that 



7! 



(po) 



2/xa 2 \ _ 2a^p 2 + q 2 



(41) 



(42) 



Equation (|4^jl is consistent with the purely classical and geometrical result from Eq. (JHHJ) - 

We can visualize the appearance of these closed orbits in the time development of the 
quantum solutions of the Schrodinger equation by constructing 2D Gaussian wave packets 
of the form 

ip G (x,y;t = 0) = ip G (x;xoiPoxib)ip G (y;yo,Poyib), (43) 



where 



ip G (x;x ,Po x ,b) 



1 



e ~(x-x ) 2 /2b 2 e ipox(x-x )/H 



(44) 



A similar expression holds for 4> G (y,yo,Po y ,b). Once again, as long as the initial location, 
(xq, yo), is far away from the edges of the potential well, such a Gaussian form can be easily 
expanded in terms of the eigenstates. The various position and momentum space expectation 
values are given by analogs of Eqs. (J1J) and ©, and the expectation value of total energy is 
now 

1 _ 1 . _ _ fi 2 

(45) 



(E) = i- (p 2 + p 2 ) = i- [ P0x 2 + p 0y 2 + ^] . 
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The expansion coefficients also exhibit separability with 

oo oo 

M*> V\ t = 0) = E E ^^WK^HM^^^ ( 46 ) 

(47) 



n x =l n v =l 



x — J- lt>y 
OO 

)t/fi 



= 4> G (x,t)ip G (y,t). 
with an equally simple expression for the autocorrelation function, namely 

A(t) = A x (t)A y (t). (48) 

We illustrate the time-dependence of such wave packets by plotting |A(t)| 2 = 
\A x (t)A y (t)\ 2 versus t in Fig. 4. The results shown are for wave packets character- 
ized by initial positions (x ,y ) = (a/2, a/2) and initial momenta given by (po x ,Po y ) = 
(po cos(#),po sin(0)) where po = 4007T and 9 = tan _1 (poy/pox); we also use Ax$ = Ay = 0.05 
and the same physical parameters as in Eq. (J25|) . With these values, the classical period (for 
the simplest back-and-forth motion), T cl = 2a/ v , the spreading time, t (from Eq. (fT2|0. and 
the revival time, T rcv (from Eq. (jHHO, are given by T c \ = 2fxa/p = l/(4007r) ~ 0.8 x 10 -3 , 
t = (2fi/h)Ax 2 = (0.05) 2 = 2.5 x 10" 3 , and T rev = 4^ia 2 /(hn) = 2/n w 0.64. The wave 
packet will then exhibit a reasonable number of classical periods before significant spreading 
occurs, with the revival time scale being much larger than both (which is a typical hierarchy 
for the time scales for physical revivals, such as seen in Rydberg atom states.—) 

We show |v4(t)| 2 versus t in Fig. 4 for increasing values of 9 up to 45° (about which the 
pattern of classical orbits and autocorrelation function data are symmetric). We indicate 
the expected location of the classical closed or periodic orbits by stars at locations of values 
of T cl p /t = \Jp 2 + q 2 and 9 = tan _1 (g/p) corresponding to the values in Table I. (Images 
of some of the simplest classical closed orbits are also included at the far right for the ap- 
propriate values of 9.) We note that the classical periodicities for closed orbits are obviously 
reproduced in the time-development of the appropriate quantum wave packets. We have also 
confirmed that varying the initial location of the wave packet has no affect on the results 
shown in Fig. 4, consistent with the classical result that the existence of closed orbits depend 
only on p, q through Po y /po x = tan(#) = q/p. 
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C. Special triangular billiards 



x,y) = 


1 


[u {n) (x)u {m) (y) - u {m) (x)u {n) (y)] 


(m 7^ n) 


(49a) 


x,y) = 


i 
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[u( n )(x)u( m )(y) + w (m )(x)it( n )(2/)] 


(m ^ n) 


(49b) 


x,y) = 




»«(n)(2/), 




(49c) 



The energy eigenvalues and wavefunctions for a special 2D triangular billiard can be easily 
derived from those of the infinite square well solutions in Eq. (J35|) . The standard results for 
the square well, E(n x ,n y ) = (h 2 n 2 /2[ia 2 )(n 2 x + n 2 ) and w {nx ^ y) (x, y) = u M (x)u( ny) (y), hold 
for any integral n x ,n y > 1; for n x = n y there is a single state, while for n x ^ n y , there is a 
two-fold degeneracy. Linear combinations of these solutions can be written in the form 

W (n!m) { 

w {+) I 

W (n,n) { 

which have the same energy degeneracy, but exhibit different patterns of nodal lines. These 
alternative forms are useful because they allow one to discuss the energy eigenvalues and 
eigenfunctions of the 45°-45°-90° triangle billiar d 39 i 40 i 41 formed by "folding" the square along 
a diagonal, because the w^ n ^ (x, y) satisfy the appropriate boundary condition along the new 
hypotenuse. (We can easily see from Eq. ()49a|) that wi' m Jx, y = x) = by construction.) 
Additional foldings along diagonals are also possible,— and the energy spectrum can be 
used to analyze these CclSCS clS well. The allowed eigenvalues for this case are still given by 
E{n xi n y ) = (h 2 7i 2 /2fia 2 )(n 2 +n 2 ), which corresponds to 

h 2 Ti 2 

E{n,m) = —{n 2 + m 2 ), (50) 

but now with only a single (n x , n y ) — > (n, m) state allowed, with corresponding wavefunctions 
given by Eq. ()49aj) . but multiplied by v^2 to account for the different normalization needed 
in the smaller area billiard. 

The revival time in the 45°-45°-90° triangular billiard is still given by Eq. (JHUj) and 
localized wave packets can also be constructed using the appropriately normalized analogs 
of the wavefunctions in Eq. (}4T?|) . again provided they are kept away from any of the infinite 
wall boundaries. 

The structure of the classical closed or periodic orbits in this case is the same as for the 
square billiard, because all of the standard (p, q) orbits in the square well are simply reflected 
off the new diagonal wall (along the hypotenuse), giving rise to the same allowed orbits as 
in Eq. (JHTj) and Table I. The only new feature in the semi-classical propagation of such 
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wave packet solutions in this folded geometry 4 ^ is the existence of a special isolated closed 
orbit at 135° (one that bisects the 90° right angle, bouncing normally off the hypotenuse, 
and returning to the corner), which has path lengths that are multiples of (\/l 2 + l 2 )a/2 = 
y / 2a/2, namely half that of the standard (p,q) = (1, 1) features. When we construct wave 
packets using parameters appropriate to this geometry, we see twice as many features in 
the A(t) plot for this case because of this special classical closed orbit, but we otherwise 
reproduce the results shown in Fig. 4. 



V. EQUILATERAL TRIANGLE BILLIARDS 

It is perhaps under appreciated that the energy eigenvalues and position space wave 
functions for a particle of mass fi in an equilateral triangular (60°-60°-60°) infinite well (or 
billiard) of side a are available in closed form. They have been derived for a variety of dif- 
ferent contexts by at least four groups 42 ' 4 3i44 i 4 ^ using complementary methods of derivation. 

For definiteness in what follows, we will assume such a triangular billiard with vertices 
located at (0,0), (a/2, >/3a/2), and (—a/2, \/3a/2). The resulting energy spectrum is given 
by 

L{m } n) = -{ — J [m + n — mn) (51) 

for integral values of m and n, with the restriction that m > 2n. (In what follows, we will 
use the notation of Ref. [43, except for a minor relabeling of the energies and wavef unctions: 
we will continue to use n for the particle mass to avoid confusion with various quantum 
numbers.) For m > 2n, there are two degenerate states with different symmetry properties^ 
that can be written in the forms 



w <J.n)( x >y) 
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(m,n) v ' " / \ 9 



27r(2m — n)x\ ( 2iiny\ 
sm sm — -=— 



3\/3 L V 3a J V V3a ) 
I 'luiln — m)x\ ( 27rmv\ 

'27r(m + n)x\ f 2n (m — n)y 
sm I — — sm ' 



3a J V v 7 ^ 
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and 



w 



(+) 

(m,n) 
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2ii(2m — n)x\ f 2irny\ 
cos I sin 



3a J \ y/3a J 

'2n{2n — m)x\ ( 2Tcmy\ 

cos I sm — = — 

3a J V V^a ) 

'2n(m + n)x\ (2i{{m — n)v 
+ cos I — K — — sin ' 



(53) 



3a J \ \/3a 

We can confirm by direct differentiation that they satisfy the 2D Schrodinger equation, with 
the energy eigenvalues in Eq. (|51|) . as well as the appropriate boundary conditions. We have 
here also included the correct normalizations, because we will eventually expand Gaussian 
wave packets in these eigenstates. 

For the special case of m = 2n, there is a single non-degenerate state for each value of n 
given by 



w 



(0) 

(2n,n) 



2imx\ ( 2iTny\ 

2 cos I ) sm — -=— 

\V3a J 



. ( 47rny 



(54) 



Clearly these states satisfy 



w< $,n)(- x >y) 



s(x, y) 

(m,n) V ' y / 



+w 



(0) 

(m,n) 



(55a) 
(55b) 



and the w^ =2n n ^ (x, y) states also satisfy 



w 



w 



(+) 

(m=2n,n) 
(") 

(m=2n,n) 



[x,y) = y/2- 



w 



(0) 

(2n,n) 



[x,y) 



o. 



(56a) 
(56b) 



The pattern of energy level degeneracies and wavefunction symmetries is thus very similar 
to that for the square billiard, especially when the solutions for that problem are written in 
the form of Eq. flUJ) . 

We note that the corresponding momentum space wavefunctions, f(^ n \(Px,Py), can also 
be obtained in closed form, because the required 2D Fourier transform over the triangular 
region can be done using standard integrals involving trigonometric functions and powers, 
making this problem very useful for still another reason. 

We turn now to the time-dependence of wave packets in this geometry and consider the 
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long-term, revival time scales in this two quantum number system: 

rev 

rp(n) 

rev 
rp(m,n) 



\d 2 E/dm 2 \/2 

2nh 
\d 2 E/dn 2 \/2 

2irh 
\d 2 E/dmdn \ ' 



which gives 



rp(m) r pij l ) rp(m,n) _ rp 



9/ia 2 
Arm 



(57a) 
(57b) 
(57c) 

(58) 



Exact revivals are present in this system (just as for the square well), with a single revival 
time guaranteed for any possible initial wave packets. Thus, both the N = 3 and N = 4 
polygonal billiards (the equilateral triangle and square) exhibit similar and simple energy 
eigenvalues and exact quantum revivals. We note that special, shorter-time scale revivals 
are also possible in the equilateral triangle case for zero-momentum wave packets initially 
centered at symmetric locations within the triangular billiard,^ such as at the geometrical 
center and half-way down a bisector. 

The short-time classical periodicity of quantum wave packets in this geometry can also 
be determined from calculations of 



(59a) 
(59b) 



rp(m) _ 
1 cl — 


2tvH 


"9/ia 2 


1 


\dE/dm\ 


4fm 


(2m — n) 


rp(n) _ 
1 cl — 


2nh 


9/ia 2 


1 


\dE/dn\ 


Ahir 


(2n — m) 



The condition leading to closed orbits can then be written as 



(2m — n) T c 
(27T 



(n) 



m) 



T, 



(m) 



cl 



P 

Q 



or 
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m 



2p + q 



(60a) 
(60b) 



2q + p / 

If we substitute this condition into the energy spectrum in Eq. (JBTj) and equate the quantum 
energies with the classical kinetic energy, [IVq/2, we are led to the association of \iv\j2 with 



E(m, n) and 



E(m, n) 



16tt 5 



h 2 



2/ia s 



2 2 

m + m 



2p + q 
2q + p 



m 



2p + q 
2q + p 



(61) 
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or 



UttS j ~ L (2q + P ) 2 J ' 1 j 



This identification implies that 



m = ( ^ ) - y (2g + p) (63a) 
V 4:7rh J y^Jp 2 +pq + q 2 



- - ( 2/a T) . (2p+q) ( 63b ) 

V 4vrfi, J ^ZJj? + pq + q 2 



The period for classical, closed/periodic orbits is then given by 



rp(po) ™(m) aV3^/p 2 +pq + q 2 L(p,q) 

T « =PT - = V = (64) 

where 

L(p, q) = d{p, q) = aV3\/p 2 + pq + q 2 (65) 

are the corresponding path lengths for periodic orbits. The possible classical closed or 
periodic orbits can be derived 32 from a geometrical construction (involving tiling of the 2D 
plane, just as for the square case2L2&) giving the same result. We note that the special 
cases of 2m = n and 2n = m correspond to q = and p = respectively, both of which 
give L(p, q) = y/3a. In these cases, one of the classical periods from Eq. (}59|) formally goes 
to infinity, which can be understood classically from the corresponding path length, which 
corresponds to periodic back and forth motion from one corner, along a bisector, to the 
opposite side, but with no repetition in the complementary direction. (See Ref. |43j for a 
derivation of the quantized energies from which this effect also can be inferred.) 

Because of the (relatively) simple form of the allowed wavefunctions in Eqs. ()52j) and 
(153)1 . we can evaluate the expansion coefficients for any 2D Gaussian wave packet of the 
form in Eq. (j4*3*J) by extending the region of integration from the (finite) triangular region to 
the entire 2D space, as long as the initial wave packet is far away from any of the walls.— 
The required Gaussian-type integrals have the forms 



-co 



<■'■ < - (0S(— )dx= e iCx /a e -b^C/a+ P0 /hy/2 (gg) 

& 2 - 

iCx /a-b 2 (-C/a+p /h) 2 /2 



J,po{x-x ) / h ^-(x-x ) 2 /2b 2 

oo 



and 



-co 



Jp (x-x )/h -(x-x ) 2 /2b 2 



,Cx^ , bV2n 



i -( - sm(— )dx = -^r- e iC*o/a e -b 2 (C/a +Po /hy/2 (g ? ) 

_ e -iCx /a e -b 2 (-C/a+p /h) 2 /2 
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with similar expressions for the y-integrations. The resulting closed form expressions for the 
expansion coefficients can be used in calculations of the autocorrelation function to illustrate 
the long-time revival structure of the wave packets, as well as the short-time, semi-classical 
propagation giving rise to closed orbits of the form in Eq. (J64j) ; the analogs of Figs. 2 and 4 
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for the equilateral triangle billiard have been presented in Ref. 

A folding of the equilateral (60°-60°-60°) triangle along a bisector yields another special 
triangular geometry, namely a 30°-60°-90° right triangle. The energy eigenvalues and eigen- 
functions for this case can also be trivially obtained from those of Eq. ()52|) as they satisfy the 
new boundary condition along the fold. Such a system will have the same energy eigenvalues 
as in Eq. (J51j) (but with only one possible (m, n) combination with m > 2n) and the same 
common revival time, T rev , as in Eq. (J5S|) . Wave packet solutions can also be constructed (re- 
membering to include an additional factor of to account for the normalization difference) 
( 

(m.n) v 



from the w\ n \(x, y) in Eq. (J52|l . 



VI. CIRCULAR BILLIARDS 

A. Eigenstate solutions and connections to classical closed orbits 

We finally turn our attention to the iV — > oo limit of the iV-sided regular polygonal 
billiard, namely the circular infinite well. (We note that the problem of scattering from an 



46l ) The potential 



infinite cylindrical barrier in two dimensions as been considered in Ref. 
for this problem can be defined by 

{0 for r < R 
(68) 
oo for r > R 

The (unnormalized) solutions of the corresponding 2D Schrodinger equation are given by 

u {m) {r,B) = J lml (kr)e imt> , (69) 

where the quantized angular momentum values are given by L z = mh for m = 0, ±1, ±2, . . . 
and the J\ m \(kr) are the (regular) Bessel functions of order \m\. 

The wavenumber, k, is related to the energy via k = ^2^E/h 2 and the energy eigenvalues 
are quantized by the boundary conditions at the infinite wall at r — R, namely J\ m \(z = 
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kR) = 0. The quantized energies are then given by 

E M = ^ R2 , (70) 

where zr m>n \ denotes the zeros of the Bessel function of order \m\, and n r counts the number 
of radial nodes. 

The energy spectrum is doubly degenerate for \m\ ^ corresponding to the equivalence 
of clockwise and counter-clockwise (m > and m < 0) motion. We therefore see a pattern 
of degeneracies very similar to that of both the square and equilateral triangle wells, with 
two equal energy states for each (|m| > 0,n r ) value, and a single one for each (m = 0,n r ). 
Because the quantum number dependence of the energy eigenvalues is the determining fac- 
tor in the structure of wave packet revivals, we need to examine the m, n r dependence of 

E(jn,n r ) ^ [^(m,n r )] • 

As a first approximation, we can look at the large z behavior of the Bessel function 
solutions 47 for fixed values of \m\, namely 

«)-^co S ( 2 -M_!) + ... (71) 
With this approximation, the zeros would be given by <2( m ,„ r ) — — f ~ (2n r + 1)|, or 

I TTli I 3 

Z(m,n r ) ~ fa + — + -)tT. (72) 

If this result were exact, the quantized energies would depend quadratically on two quantum 
numbers, and there would be exact wave packet revivals, just as for the square or equilateral 
triangle billiards. However, there are important corrections to this first-order result, which 
means that the Bessel function zeros are not given by exact integral values. So although 
there are not exact quantum revivals in the general case, approximate revivals are present, 
most especially for zero-momentum, central ((x ,y ) = (0,0)) wave packets (constructed 
from purely m = states). We refer the interested reader to Ref. |3l| for details. 

We can, however, examine the short-time dependence leading to classical closed orbits 
by making use of the WKB approximation to evaluate the quantized energies and their 
m and n r dependence. If we first quantize the angular variable to find that the angular 
momentum is given by L z = mh, we note that in the radial direction the particle moves 
freely up to the infinite wall air — R, but is subject to an effective centrifugal potential given 
by Ks( r ) = Ll/2fir 2 = (mh) 2 /2[ir 2 . A classical particle cannot penetrate this centrifugal 
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barrier and therefore has an inner radius (distance of closest approach) given by Ves(Rn 
m 2 h 2 /{2 M i4 in ) = E, or 



\m\h 

VP 



(73) 



\m\R/z and E = h 2 z 2 /(2/jR 2 ) in terms 



We can also write Eq. (|75j) in the useful form R mir 
of the dimensionless variable, z. 

The WKB quantization condition on the radial variable r is then given by 

-R 



k r (r) dr = (n r + C L + C R )n, 



Rn 



where 



kJr) = 



2fiE 



d2 

-y min 



(74) 



(75) 



h 2 V r z 

The matching coefficients^ are given by Cl = 1/4 and Cr = 1/2 which are appropriate 
for linear boundaries (at the inner centrifugal barrier) and hard or infinite wall boundaries 
(such as at r = R), respectively. (See Appendix A for more discussion of WKB methods.) 
The WKB energy quantization condition for the quantized energies in terms of n r explicitly 
and \m\ implicitly can then be written in the form 




>E-^dr 



2/ir 2 



2fj,E 



R 



R 2 ■ 
1 ^dr 



3/4K 



(76) 



which defines the quantized energies, E = £ , ( m „ r ), implicitly in terms of m,n r . To obtain 
the partial derivatives that are necessary to evaluate the classical periods in Eq. (|32|) . we 
can then differentiate this condition with respect to both quantum numbers to obtain 



2h 2 



dr 



fl mi „ JE - m 2 h 2 /2/j,r 2 



2h 2 
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dr 



fdE 



dE 

dn r 

m\h 2 
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/ Rmin yjE - m 2 h 2 /2fir 2 \dm fir 2 
The condition for periodic orbits from Eq. (|34j) can then be written as 



(77a) 
(77b) 



q 
p 



T, 



(nr) 



T, 



(m) 



cl 



\dE/dm\ ( \m\h \ 
\dE/dn r \ ~ y-K^EJ 
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(78) 



If we evaluate the integral and use i? n 



\m 



\h/y/2jlE, we find that 



Q 
P 



1 _i/ R s 
_ sec ( ) 



7T 



(79a) 



or 
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Rmin(P,q) = Rmin = # COS (79b) 

as the condition on periodic orbits. To find the classical period for such closed orbits, we 
note that 

t (p°) _ „ T K) _ r 27rft i _ ( ~ [W 2 ^ IT _ 2pRsm{nq/p) _ L{p, q) 

where we identify vq = \JlEj n with the classical speed. The condition in Eq. (J79)l and 
the resulting path lengths/periods in Eq. ()80|) can be obtained from strictly geometrical 
arguments .^ 1 ^ Just as for the square and equilateral triangle cases, information on the 
classical closed orbits is clearly encoded in the energy eigenvalue spectrum for the circular 
well using Eqs. (jH2^1 and (JBH). The values of these path lengths, L(p,q) = 2pRsm(7rq/p), 
and the corresponding distances of closest approach, R min /R = Rcos{nq/p), for many of 
the lowest-lying closed or periodic orbits in the circular well are collected in Table II. 

B. Wave packet construction in the circular billiard 

To examine the construction of Gaussian wave packets in this geometry, we note that 
any arbitrary initial wave packet, tp(r,8), in the circular billiard can be expanded in the 
normalized eigenstates of the form 

wV,„ r) (r,#) = [N (m>nr) J\ m \(k^ nr) r)](-=e ime ), (81) 

V Z7T 

where 

[N {m , nr) ] 2 r[J lml (kr)] 2 dr = l, (82) 
Jo 

with expansion coefficients given by 

a(m,n r ) = (ip(r,0;t = 0)\w {mjrir) ), (83) 

which satisfy 

+oo oo 

E J2\ a (m,nr)\ 2 = l- (84) 
m=— oo n r =0 

Because the eigenstates are no longer simple trigonometric functions, the required normal- 
ization and overlap integrals must be done numerically, making the process of wavepacket 
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construction much more involved. The expectation value of the energy in this potential well 
is given by 

+00 00 +00 00 / f% 2 ^Z ] 2 \ 

(E)= E E \ a (m,n r )\ 2 E {mi n r) = E E l°(«».«r)| 2 ( 2 uB? J ' ^ 

m=— 00 n r =0 m=— 00 n r =0 \ r" / 

and this constraint is another useful one for numerical checks. 

In this more symmetrical geometry, we also can evaluate expectation values of powers of 
the angular momentum, 1^, which is the other important conserved quantity. We find that 

+00 00 

(4 fe >= E EKwv)| 2 M^- (86) 

m=— 00 n r =0 

The subsequent time-dependence of any such wave packet is then given by 



-00 00 



#",0;t) = E Y a (™^)nrn,n r )(r,0)e- iE ^/ n , (87) 



m=— 00 n r =0 

and the autocorrelation function is given by^2, 



-00 00 



A(t) = (iKr,e;t)\if>{r,0,O))= E E l a KMl 2 ^ K " r)!M • (88) 

m=— oo n r =0 

We now focus on the specific Gaussian form in Eq. (|43|). besides the expectation values 
calculated above for such a form. We can explicitly evaluate the conserved powers of the 
angular momentum, namely 

(L z ) = (xp y - yp x ) = (x)(p y ) - (y)(p x ) = x p 0y - y p 0x , (89) 

and 

b 2 h 2 
(L 2 Z ) = (x oPoy - y p 0x ) 2 + - [(p 0x ) 2 + (p 0y ) 2 ] + — [(x ) 2 + (y ) 2 ] , (90) 

so that the spread (or uncertainty) in the angular momentum is given by 



b 2 h 2 
(Am)h = AL Z = yj- [(p 0x ) 2 + (p 0y ) 2 } + — [(x ) 2 + (y ) 2 }. (91) 

This form can be understood intuitively from a simple error propagation argument that 
would give the error in the product L ~ rp to be (AL/L) 2 = (Ar/r) 2 + (Ap/p) 2 or 

AL = vV(Ar) 2 + r 2 {Ap) 2 . (92) 
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We note that even for wavepackets for which the expectation value of the angular momen- 
tum vanishes, it may still be necessary to include m/fl components; for example, a wave 
packet with (x ,y ) = (0,0) and (po x ,Poy) — (0,Po)> will have a necessary spread in angular 
momentum values given by Am = bpo/\/2h, which increases linearly with the initial momen- 
tum. For such cases, because (L z ) = 0, we must have relations such as |a( m)Jlr )| = \a(- m ,n r )\ 
to ensure that Eq. (|8Hj) is satisfied. 

Because of the required numerical evaluation of the expansion coefficients, it is useful to be 
able to compare the general results for (E) and (Li ) in Eqs. (|55)l and (JHtjj) with the specific 
results for the Gaussian in Eqs. 1)45)) . ()89j) . and (|9()j) . As an example of such an expansion, we 
have numerically evaluated the at m>nr ) for Gaussian wave packets characterized by yo, Pq x = 0, 
Po y = 100, and xq/R = 0.0, —0.5 and 0.70, that is, with the same total (kinetic) energy, 
but different values of the angular momentum; the same physical parameters of Eq. (|23|) 
are used, along with R = 1. The various expansion probabilities, |a( m ,n r )| 2 , are sorted by 
magnitude, and then summed to confirm that the numerical evaluation has captured all of 
the probability density. We can then illustrate which regions of the (m, n r ) space are most 
populated in the construction of each wave packet and the results are shown in Fig. 5. In 
each case, the (approximately elliptical) inner shaded areas correspond to regions in (m, n r ) 
space containing 68% (inner dark region) of the total probability, while the outer dark 
regions surround 99.7% of the ^(mrv) \ a (™.,n r )\ 2 ■ We note that the required spread in m 
values increases for those states with increasing values of | (L) \ = xopo y . 

We can also follow the short-time development of Gaussian wave packets in the circu- 
lar billiard for various initial conditions to look for evidence of the classical periodicities 
described by Eqs. (|T9"j) and (JHUj) . In this case we construct packets with y Q = 0, p 0x = 0, 
Po y = 100 and various values of Xq = R m i n - (We use a much smaller value of po than in 
earlier calculations for the square well because we are restricted to evaluating the expansion 
coefficients numerically which becomes prohibitively time-consuming for larger values: the 
repeated evaluation of high-order Bessel functions is numerically intensive in the programs 
we use and is the limiting factor. This unavoidable restriction means that the spreading 
time, to, is of the same order as the classical periods shown.) We plot in Fig. 6 the autocor- 
relation function |^4(t)| 2 versus t/r, where r = R/vq. We note that isolated features in A(t) 
are present at the values of Jr = 2psm(irq/p) and R min / R = cos(irq/p) values given in 
Table II, corresponding to classical closed orbits in this geometry. 
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We can extend this analysis to the half-circle billiard, obtained by folding across a diam- 
eter in the same way as we have done for the half-square and half-equilateral triangle. In 
this case, we use a single copy of the m > eigenstates.— 

VII. CONCLUDING REMARKS 

We have examined the quantum mechanical time-development of localized (Gaussian- 
like) wave packets in three billiard systems, with square, equilateral triangle, and circular 
footprints. We have used the well-known (and not-so-well-known) energy eigenvalue spec- 
trum in each case to discuss the classical periodicity and quantum revival time scales in such 
systems using the autocorrelation function as a probe. In the first two cases, we used the cor- 
responding energy eigenfunctions to calculate closed form, analytical approximations for the 
expansion coefficients of the Gaussian wave packets, suitable for fast numerical evaluations 
for arbitrarily large momentum values, provided the initial wave packets are well-localized 
within the well and reasonably far from any infinite wall boundary. For the circular well, 
we have performed a similar analysis, using a WKB analysis to demonstrate the expected 
short-time classical periodicities and a numerical approach to evaluate A(t). 

In each case, we have also been able to consider half-well geometries, obtained by folding 
of the original footprint along an obvious axis of symmetry for which solutions were available 
as a subset of those for the full-well. 

Further extensions of these studies could include more detailed examinations of rectan- 
gular geometries, or even three-dimensional (3D) parallelepipeds. The case of an annular 
billiard (with both outer and inner infinite walls) 4 " can be handled in much the same way 
as described here for the 2D circular well, and the same WKB approach can be used to 
extract the pattern of allowed classical closed orbits, which is much richer than in the simple 
circular billiard. 

The spherical billiard also can be discussed with many of the same techniques. The 
energy eigenstates now consist of products of spherical harmonics, Yi tm (9,4>), instead of 
the simple e imd angular states, with spherical Bessel functions for the radial component. 
We might expect that the resulting change in the angular momentum term from two to 
three dimensions would be given by the substitution m 2 — > /(/ + 1), but it turns out that 
for a WKB analysis it is more appropriate 48 to use (I + 1/2) 2 instead of 1(1 + 1) (the so- 
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called Langer modification.—) This identification still gives the same result for the classical 
periods, which is not surprising because angular momentum conservation implies that all the 
allowed 3D classical orbits will be planar and reduce to the 2D case. This effect can be seen 
rather generally in 3D central potential systems where we know that the quantized energy 
eigenvalues do not depend on the m quantum number, namely the E(n r ,l,m) = E(n r ,l) 
are functions of n r and I alone. This independence implies that the classical periodicity 
corresponding to the m quantum number in Eq. (}3*2*|) will be = 2-nh/\dE / dm\ — ► oo 
which becomes irrelevant, leaving only and to beat against each other to provide 
closed orbits. 

In the context of the spherical billiard, we note that the spherical Bessel functions, ji(z), 
can be related to the ones considered here for the 2D case via ji{z) = \Jit l*lzJi+\ii{z). 
For purely central (no angular momentum, s-wave) Gaussian wave packets (therefore with 
initial position at (x , yo, z ) = (0,0,0) and with vanishing momentum), the required / = 
spherical Bessel function is proportional to sin(z)/z, which does have exact integral values 
of z — kR = nir. For this case only, there will be exact quantum revivals in the spherical 
billiard. 

The solutions for the equilateral triangle case (N = 3 regular polygon) form a subset 
of those required for a complete study of the (N = 6) hexagonal billiard and that case is 
currently under investigation. It is an interesting and open question as to whether there 
are other N values (besides N = 3 and 4) for which there exist exact quantum revivals 
as well as well as how the N — > oo limit of the circular billiard (where only approximate 
revivals are present) is reached. Approximation methods appropriate for such quantum 
billiard systems^ may be useful in answering such questions. 

We also can extend the notions explored here dealing with classical periodicity to non- 
billiard systems. For example, the extension of Eqs. (|32|) and (|34|) to three dimensions is 
obvious. As mentioned, for central potentials the classical period is irrelevant, so that 
closed or periodic orbits are the result of the condition pT^ r ' = = qT® for the (n r , I) 
quantum numbers. It is a familiar result of classical mechanics (Bertrand's theore m^^ 4 ) 
that for the 3D power law potentials (V(r) oc r k ), only the harmonic oscillator (k = 2) and 
Coulomb potentials (k = —1) give rise to closed orbits for all bound state trajectories.— It 
is straightforward to exhibit a similar statement for the quantum versions of these systems 
in terms of their expected classical periods, T^°\ 
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APPENDIX A: WKB ANALYSIS OF CLASSICAL PERIODS 

As a comparison with the formalism leading to Eq. (|15j) . an alternative derivation of the 
expression for the classical periodicity of a bound state system in terms of its quantized 
energy levels can be obtained from a simple WKB argument. For a particle of fixed energy 
E in a bound state potential, V(x), we have E = fiv(x) 2 /2 + V(x) and the short time, dt, 
required to traverse a distance dx can be obtained from this energy conservation connection 
and integrated over a single cycle to obtain the classical period via 

Jl_te^ (Ala) 
\2y/E-V(xY 
or 

T = 2 ^ dt = 2 X /| I" - dx (Alb) 
The WKB quantization condition for this potential can be written in the form^ 

nb 

% ^E n - V(x) dx = (n + C L + C R )nh (A2) 

J a 

in terms of the matching coefficients Cl and Cr. Recall that Cl,r = 1/4 at linear or 
smooth turning points where the functions are matched smoothly onto Airy solutions, while 
Cl,r — 1/2 at infinite wall type boundaries where the wavefunction must vanish,- This 
difference is sometimes described as being due to the fact that WKB wavefunctions can 
penetrate roughly one-eighth of a wavelength into the classically disallowed region, provided 
the barrier is not infinitely high.— 

Equation (jA2|) can be differentiated implicitly with respect to the quantum number n to 
obtain 

f— f b \dE n /dn\ dx 

Equation (jA3j) can be related to the classical period in Eq. (jAlbjl to give 

r n = a/27 ["(En ~ V(x)Y^ dx = (A4) 
J a \dE n /dn\ 
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as in Eq. (|15|). The most obvious example of such a connection is for the harmonic oscillator, 
where the WKB condition gives the exact eigenvalues, E n — (n+ 1/2) fku, and the classical 
period from Eq. (jA4|) is 

2nh _ 2nh _ 2n 
Tn ~ \dE/dn\ ~~h^~~ { ' 

as expected. 



APPENDIX B: EXPANSION COEFFICIENTS FOR THE ID INFINITE WELL 
IN MOMENTUM SPACE 

Although it is typical to discuss the expansion coefficients for a Gaussian wave packet 
in the ID infinite well in position space, as in Eqs. (J7|) or (J26)) . it is instructive to examine 
the same problem in momentum space. As mentioned in Sec. 1111} the requirement of the 
necessary matching of the oscillatory e~ tPxX ^ h factor in ^(x, 0) with the corresponding 
e ±mnx/a f ac t ors f r0 m the bound state u n (x) is contained explicitly in the exp(— b 2 (po ± 
rnrh/a) 2 /2h 2 ) terms in Eq. (f2*U|) . There must also be an appropriate overlap between the 
envelope of iPg{x, 0) contained in the exp(— (x — Xq) 2 /2b 2 ) factor and the extent of the bound 
state wavefunctions. That is, it must be inside the well, as illustrated in Fig. 1. 

The form of the momentum- matching can be qualitatively understood in the following 
way. The oscillatory parts of the integrals required for the evaluation of a n will contain 
factors such as exp(i(p ± mrh/a)x/h) to be integrated over x. An integral over all space 
would give rise to a ^-function term of the form 5(\po±mih/a}/h), because of the cancellations 
arising from the rapid oscillations of the integrand. However, the integrals are effectively cut 
off by the spatial extent of the initial wave packet, a width of order Ax ~ b, which instead 
gives a sharply peaked function of \po ± mih/a\b/h, as seen in Eq. (j2El)- 

The same information must be encoded in a momentum space representation of the 
evaluation of the a n expansion coefficients. The required Fourier transform 

MP) = -i= r u n {x) e-^ h dx (Bl) 



o 



V2nh 

of the bound-state wavefunction u n (x) = -^/■|sin( n ^) is given by 
a 



V 'Ait ha 
at 



V 'An ha 



' ^i(n-K—pa/ h) g—i(nir+pa/h) 

(nn — pa/h) ^ (mr + pa/h) 

c inn Sin K n71 - WW 2 ] _ e -in7r Sin [( n7r + WO/2] 



(B2) 



[(mr - pa/h) /2] [(mr + pa/h) /2] 

29 



—ipa/2h 



while the initial Gaussian from Eq. ((TJ is 



<j> G {p) = < H= e- ai{ P- pa ^/ 2 e - ipxo/h . (B3) 
V V n 

The (exact) expansion coefficients now have the form 

X) 

K(p)MP,0)dp (B4) 



The matching in momentum space now comes from the requirement that the peak in <Pg{p) 
(near po) must match with those in <f) n (p) (near p = ±nirh/a). The information on whether 
the initial wave packet is contained inside the (0, a) infinite well now arises from the phase 
information contained in the oscillatory parts of the momentum space versions: specifically, 
the a n integral in Eq. (|B4|) contains terms such as 



[<f>*n(p)][<t>G(p, 0)] OC [e"W2?i]*[ e -^oA] = e ip/h(a/2-x )^ ^ 

which, if integrated by itself over all p-space, would give rise to a S(xq — a/2) term. However, 
this singular behavior is again softened to be a sharply peaked function of \xq — a/2\, due 
to the finite extent of the p-integrals. It should not be surprising that the roles played 
by the envelope/overlap requirement versus that played by the oscillatory information are 
complementary in the x- versus p-space descriptions of the integrals in Eqs. (J21j) an d (jB4|) 
We can easily generalize this last result to the case of an infinite well of width a, but 
located at an arbitrary location along the ID axis defined over the range (d, d+ a). In that 
case, the energies are unchanged, and the position space eigenfunctions are simply shifted 
to 

u n (x) — u n (x — d) — \ - sin(— — — — ). (B6) 

V o a 

The Fourier transform to obtain the corresponding momentum space eigenstates goes 
through as before with a simple change of integration variables, x — d —>■ y, giving 

Up) = Mp)e- ipd/ \ (B7) 

in which case the appropriate phase factor, analogous to that in Eq. (jB5|) . is exp(ip/h(d + 
a/2—xo)); values of Xq beyond the boundaries of the new well lead to exponential suppression. 

APPENDIX C: PROBLEMS 

We offer some suggested problems that are motivated by the suggestions in Sec. IVHI 
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Problem 1. Extend the results of Sec. HVBI for a square to a rectangular billiard with sides 



L x x L y . Calculate the allowed energies and the c 



exact and/or fractional revivals. (See Refs. 



35 and 



assical periods. Discuss the patterns of 



36 for details. ' 



Problem 2. Extend the relationship in Eq. ()34|) to a system with three quantum numbers 
to find the allowed classical periods. Apply this extension of Eq. to the system of a 
cubical box of side L, generalizing the result of Eq. (}4*2j) . 

Problem 3. Visualize some of the closed orbits for the equilateral triangle geometry. (See 
Ref. I32J for examples.) 

Problem 4. Generalize the WKB results of Sec. IVI Al to show how the closed orbits (periods 
and path lengths) in an annular infinite well are obtained. How is Eq. (J76|) changed, including 
the limits of integration and the matching coefficients, Cx and Cjp. (See Ref. nd for details.) 

Problem 5. Plot the equivalent of Fig. 5 for Gaussian wave packets in a square well, initially 
centered at the origin, with initial momentum components {pox,Poy) = (po cos(9),po sin(0)) 
for various values of 9. That is, evaluate the expansion coefficients, d(n x ,n y ) = a n x o>n y , and 
the corresponding probabilities, and plot where the probability is distributed in the (n x ,n y ) 
plane. 
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TABLES 



9 (deg) 


period J c j /r 


(p,q) 


(deg) 


period i cl /t 


(p,q) 


0.00 


i r\(~\ o on o oo 

1.00, 2.00, 3.00, . . . 


( 1 0\ 

(1,0) 


23.96 


o o r 

9.85 


/O /I \ 

(9, 4) 


6.34 


9.06 


/o 1 \ 

( 9 > 1) 


26.57 


O O A A /IT d 71 O O A 

2.2% 4.47, 6.71, 8.94 


/O 1 \ 

(2, 1) 


7 1 O 

7.13 


8.06 


/" O 1 \ 

( 8 , !) 


OO 'T /I 

29.74 


o or* 

8.06 


(5, 6) 


8.13 


7.07 


/ T 1 \ 

(7, !) 


on n/ 1 

30.96 


c o o 

5.83 




o i 

9.46 


f* O o 

6.08 


(6, 1) 


OO oo 

32.00 


O /I O 

9.43 


(8, 5) 


11.31 


r 1 o 

5.10 


/r i \ 

( 5 > 1) 


o o r* o 

33.69 


3.61, 7.21 


/ O O \ 

(3,2) 


1 o ro 

12.53 


o oo 

9.22 


( 9 , 2 ) 


or r /i 

35.54 


8.60 


( 7 , 5 ) 


1 A O A 

14.04 


4.12, 8.25 


(4,1) 


36.87 


r oo 1 o oo 

5.00, 10.00 


/ /I o \ 

(4,3) 


15.95 


I.Zo 


/ 7 OA 

(7,2) 


3o.oo 


0.40 


( 5 ,4) 


18.43 


3.16,6.32,9.49 


(3,1) 


39.81 


7.81 


(6,5) 


20.56 


8.54 


(8,3) 


40.60 


9.22 


(7,6) 


21.80 


5.38 


(5,2) 


45.00 


1.41,2.83,4.24,... 


(1,1) 


23.20 


7.62 


(7,3) 









Table I: The periods for the classical closed orbits, T cl /r = a/p 2 + q 2 , where r = 2a/ vq, 
and the corresponding initial angles, tan(0) = g/p, for the square billiard. All periodic orbits 
with T c ( f o) /r < 10 are included. Values for 45° < 9 < 90° are the same as those for 90° — 9. 
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Tj iR = It) sin ( 7T(7 / 7) I 


/ hi 

cos^Trq jp) 




TjlR = Id sin ( Tin vt ) 


hi / hi 

cos^irq /p) 


(2,1) 


4.00,8.00, 12.00, 16.00, . . . 


0.00 


(6,3) 


12.00, . . . 


0.00 


(3,1) 


5.20, 10.39, .... 


0.50 


(7,3) 


13.65, . . . 


0.22 


(4,1) 


5.66, 11.31, .... 


0.71 


(8,3) 


14.78, . . . 


0.38 








(9 3) 


15.59, 


0.50 


(oo 1) 


6 28 12 57 


1.00 


(10 3) 


16.18, 


0.59 








(11,3) 


16.63, . . . 


0.66 


(4,2) 


8.00, 16.00, . . . 


0.00 


(12,3) 


16.97, . . . 


0.71 


(5,2) 


9.51, . . . 


0.31 


(13,3) 


17.24, . . . 


0.75 


(6,2) 


10.39, . . . 


0.50 








(7,2) 


10.95,... 


0.62 


(oo,3) 


18.85,... 


1.00 














(oo,2) 


12.57,... 


1.00 









Table II: Tabulated values of the path length, L/R = 2psm(nq/p), from Eq. (|8T)|) and 
the distance of closest approach, R m[n / R = cos(7rg/p), from Eq. (f?9*|) for low-lying closed 
or periodic orbits for the circular well characterized by (p, q) . All closed orbits with path 
lengths satisfying L(p, q) < 20R are included, as are their low-lying recurrences. 
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Figure Captions 



Fig. 1. Plots of the expansion coefficients, Y?m=i \ a n\ 2 (bottom), and the average energy, 
'Y^ = \E n \a n \ 2 (middle), for zero-momentum Gaussian wave packets in the ID infi- 
nite well, as the central position, xq, is varied from the center (x^/a = 1/2) to beyond 
the edge of the well (x /a > 1). The expression in Eq. is used and contributions 
from the 40 lowest-lying states are included. The solid (dashed) curves correspond 
to Axo = 0.05 (0.1). As expected, the normalization decreases uniformly from unity 
as more of the initial wave packet is outside the well. The average energy, however, 
increases as the initial wave packet is placed near the edge of the well, as the expansion 
in eigenstates attempts to reproduce the sharp discontinuity at the wall (as shown for 
the Xq = 0.95 packet at the top.) 

Fig. 2. Plots of the autocorrelation function, |v4(t)| 2 versus t, over one revival time for zero- 
momentum (p = 0) Gaussian wave packets in the ID infinite well. Plots for different 
initial position values, xo, from the center (xo/a = 0.5) to near one edge (xo/a = 0.8) 
are shown. Exact revivals at T rev = Afia 2 /hn are obvious for all initial positions, but 
special exact revivals at shorter times due to obvious symmetries in the eigenstate 
expansion at Xq = a/2 (for multiples of T rev /8, one of which is highlighted in the 
rectangular area) and at x = 2a/3 (for multiples of T rev /3, highlighted in the two 
elliptical areas) are also apparent. We note that there are partial revivals for x$/a = 0.8 
(top line) for times given by 0.4T rev and 0.6T rcv (shown by the diamonds.) 

Fig. 3. Plots of the autocorrelation function, |^4(t)| 2 versus t (left column), over one revival 
time for xq = a/2 wave packets with increasing values of average momentum, p Q . For 
the bottom case (d), the classical periodicity is so short that the individual periods 
cannot be resolved, in contrast to case (c) where the periods are clearly visible. Lo- 
cations of fractional (partial) revivals are indicated by arrows for case (d) where they 
are most obvious. The corresponding values of the expansion coefficients, \a n \ 2 versus 
n, are shown on the right. 

Fig. 4. Plots of the autocorrelation function, |A(t)| 2 versus t, for the square billiard. The plots 
are over a time period to ten classical back-and-forth periods, lOr, where r = 2L/vq. 
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Plots for different values of the initial angle, tan(0) = po y /pox are shown. The stars 
indicate the positions of classical closed orbits (and recurrences) as shown in Table I. 

Fig. 5. Values of the (numerically evaluated) zeros of the Bessel function, J\ m \(z = ka) as the 
quantum number, m, which determine the energy eigenvalues for the circular well, is 
varied. The (roughly elliptical) shaded regions correspond to regions in (m, z = fc m , nr ) 
space inside which 68% (inner) and 99.7% of the J2(mn r ) \ a (m,n r )\ 2 is contained. The 
three cases shown correspond to yo,Po x = 0, p 0y = 100 and x = —0.5, 0.0, and 0.7 for 
Gaussian wave packets with the physical parameters in Eq. (|25j) . The common value 
of po y = hk = 100 is shown as a horizontal dashed line. (The faint vertical banding is 
an artifact of our plotting program.) 

Fig. 6. Plots of the autocorrelation function, |v4(t)| 2 versus t, for the circular billiard. The 
plots are for xo,po y ^ with y = p 0x = as one increases Xo/R — R m - m /R, The stars 
indicate the locations of classical closed orbits with periods and Rmm/R values given 
by Eqs. (J80|) and (|79|) . respectively as shown in Table II. Examples of closed orbit 
trajectories for several simple (p, q) combinations are also shown at the corresponding 
values of -R m in/ R- 
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la I 2 vs. n 

n 





(1 1) (4 3) (3 2) (2 1) (5 2) (3 1) (4 1) (1 0) 
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